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Two  distinct  tasks  are  covered  by  this  final  report;  (1) 
reduction  of  optical  interferometer  data  obtained  during  the  196S 
SECEDE  ill  barium  release  observations  in  Alaska,  and  (2)  optical 
interferometer  instrumentation,  observations,  and  preliminary  data 
reduction  for  the  1971  SECEDE  II  barium  release  series  in  Florida. 
Analysis  of  the  SECEDE  III  data  shows  that,  even  at  late  times,  optical 
surface  brightness  of  the  barium  ior.  clouds  (as  obtained  from  filter 
photography)  can  not  be  transformed  in  ion  column  densities  without 
applying  a  spatially  varying  correction  factor,  owing  to  sub¬ 
stantial  resonance  absorption  of  the  incoming  sunlight  as  it 
penetrates  the  cloud.  To  avoid  resonance  absorption  effects  a 
larger  (4  inch  aperture)  interferometer  coupled  to  a  16"  telescope 
was  constructed  forthe  1971  SECEDE  II  series  in  Florida.  This 
interferometer  was  operated  in  a  wavelength  scanning  mode  and  in 
a  spatial  scanning  mode  (locked  on  a  weakly  absorbed  line  frequency 
to  minimize  resonance  absorption  effects). 
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Technical  Report  Summary 

Two  distinct  tasks  are  covered  by  this  final  report;  (l)  reduction 
of  optical  interferometer  data  obtained  during  1969  Secede  III  barium 
release  observations  in  Alaska,  and  (2)  optical  interferometer  instru¬ 
mentation,  observations,  and  preliminary  data  reduction  for  the  1971 
Secede  II  barium  release  series  in  Florida. 

Optical  interferometers  of  the  Fabry-,,  erot  (F-P)  type  permit 
detailed  measurements  of  the  spectral  line  profiles  resulting  when  sunlight 
is  scattered  from  the  neutral  and  ionized  barium  atoms  released  in  the 
ionosphere  as  part  of  the  Secede  study  series.  From  these  measurements 
column  densities  and  temperatures  of  the  neutral  and  ionized  barium  can 
be  deduced.  Analysis  of  the  voluminous  data  obtained  with  a  medium  sized 
(2  inch  aperture)  F-P  interferometer  observing  the  Alaskan  Secede  III 
releases  was  impeded  by  the  fact  that  the  data  format  was  not  suitable 
for  computer  reading,  plotting  and  analysis.  As  one  of  the  tasks  of  the 
contract,  interfacing  was  constructed  to  transform  the  l/4"  magnetic  tape 
records  (containing  data  in  pulse  train  and  audio  record  form)  to  computer- 
compatible  l/2"  magnetic  tape  in  a  word  form  readable  by  means  of  a  computer 
program  developed  for  the  purpose.  The  records  of  the  most  important 
release,  event  GUM,  were  transcribed  to  1/2"  magnetic  tape  with  this  inter¬ 
facing.  Analysis  of  the  Secede  III  data  shows  that,  even  at  late  times, 
optical  surface  brightness  of  the  barium  ion  clouds  (as  obtained  from 
filter  photography)  can  not  be  transformed  into  ion  column  densities 
without  applying  a  spatially  varying  correction  factor,  owing  to  substantial 
resonance  absorption  of  the  incoming  sunlight  as  it  penetrates  the  cloud. 

In  addition,  the  interferometer  data  indicate  that  the  neutral  Ba  atoms 
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remain  quite  hot  (~2000°K)  for  tens  of  seconds  after  release;  thus  thermal 
excitation  of  the  barium  to  metastable  levels  may  contribute  significantly 
to  the  initial  cloud  ionization  (supplementing  the  two- stage  photo-ionization 
process  usually  invoked). 

To  avoid  resonance  absorption  effects  in  transforming  surface 
brightness  of  the  ion  cloud  to  ion  column  densities,  a  larger  (4  inch 
aperture)  interferometer  coupled  to  a  l6"  telescope  (to  improve  spatial 
resolution)  was  constructed  for  the  1971  Secede  II  series  in  Florida. 

The  interferometer  was  operated  in  a  wavelength  scanning  mode  (to  determine 
neutral  and  ionized  barium  atom  temperatures)  and  in  a  spatial  scanning 
mode  (locked  on  a  weakly  absorbed  line  frequency)  to  minimize  resonance 
absorption  falsification  of  the  ion  column  densities  deduced  from  the  ion 
cloud  surface  brightness.  The  time  development  of  striations  in  the 
Secede  II  ion  clouds  was  detected  with  the  instrument  and  will  be  correlated 
with  rf  propagation  experiments  from  a  beacon  passing  behind  the  ion  cloud 
(whose  position  is  recorded  on  our  videotape  records  obtained  with  a  bore- 
sighted  TV  camera).  The  ion  line  profiles  provide  correction  factors  for 
filter  photography  of  the  cloud  which  permits  relation  of  their  measured 
surface  brightnesses  to  the  ion  column  densities  at  various  positions  in 
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I .  Introduction 


This  report  summarizes  our 'efforts  in  developing  apparatus  for 
optical  observations  of  barium  releases  in  the  Secede  series  and  for  optical 
detection  of  ionospheric  electron  heating  produced  by  a  high-powered  ground 
based  transmitter,  ©ur  progress  in  analysis  of  the  Secede  and  Ionospheric 
Modification  data  is  described  and  the  results  are  discussed.  We-sha-li 
first  -present  -the  analysis  of  our  Fabry-Perot  interferometer  data  from  the 
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Secede  III  barium  release  series  in  Alaska,  1969?  Then  we  shall  describe 
the  improved  Fabry-Perot  apparatus  developed  for  the  Secede  II  series  at 
Eglin,  Florida,  1971; and  the  first-look  at  some  of  the  data  obtained  with 
our  instrument.  Finally,  we  shall  present  the  results  obtained  with  our 
spatial-scanning  photometer  used  in  conjunction  with  the  Boulder  Iono¬ 
spheric  Modification  Experiment  in  1970  to  detect  electron  heating  in  the 
F  region. 

II.  i  Secede  III  Data  Handling  and  Analysis 
A.  Data  Handling 

A  45  mm  aperture  Fabry-Perot  interferometer  which  we  had  developed 
for  auroral  studies  was  pressed- into  service  on  short  notice  for  the 
Secede  III  observations  by  the  addition  of  a  reasonably  accurate  (~  1  milli- 
radian)  pointing  head,  control  and  data  processing  electron? cs,  4-channels 
of  Brush  chart  recorders  to  present  data  in  real  time,  and  a  4-track  l/4" 
magnetic  tape  recorder,  two  tracks  recording  instrumental  data  and  two 
tracks  recording  comments  of  the  optical  pointing  operator  and  of  the 
electronics  console  operator.  Figure  1  shows  a  sample  of  the  real  time 
data  obtained  from  the  Brush  records.  Channel  4  shows  the  cyclical 
pressure  scan  (corresponding  to  wavelength  scan  of  the  spectral  line), 
the  pulse  at  the  midpoint  denoting  a  filter  change  from  the  X4934  ion 
line  to  the  X5535  neutral  line  (the  two  lines  were  scanned  sequentially). 


Channel  9  is  the  D/A  output  of  the  line  profile  scanned  in  a  single  pass, 
while  Channel  2  is  the  line  profile  accumulated  in  a  multi-channel  memory 
after  a  number  of  passes.  Channel  7  is  the  integrated  line  intensity  of 
the  \4934  and  X5535  lines  obtained  in  a  monitor  channel  (one  D/A  overrun 
is  apparent  midway  in  the  record).  The  instrument  obtained  useful  data  on 
most  of  the  rel:ases;  however  reduction  of  the  real  time  records  was  very 
slow  and  laborious  (the  details  of  this  data  reduction  are  given  in 
Section  116),  and  the  l/4"  magnetic  tape  records  were  not  in  a  form 
compatible  for  computer-assisted  reduction. 

As  part  of  the  present  contract,  interfacing  was  constructed  to 
transform  the  pulse  train  data  on  the  1/4"  tape  in  2  stages,  first  to 
1/2"  magnetic  tape  in  an  intermediate  form,  and  then  to  another  1/2" 
magnetic  tape  in  a  word  form  suitable  for  reading  by  a  computer 
(instructed  by  an  appropriate  program  that  was  developed  for  the  purpose). 
This  interfacing  device  was  used  to  transcribe  the  GUM  release  inter¬ 
ferometer  data  (evidently  the  most  interesting  release  in  the  series); 
this  l/2"  magnetic  tape  record  can  be  made  available  on  request  to 
contractors  in  need  of  the  information  for  Secede  III  analyses.  However, 
transcription  of  the  single  record  took  about  one  month,  since  the  trans¬ 
criber  must  listen  to  the  audio  tracks  to  note  event  occurrences ,  generate 
file  marks,  and  prepare  a  log  to  correlate  all  of  the  data  inputs.  For 
example,  boresight  camera  pictures  were  taken  at  varying  intervals  and 
with  various  exposures.  Pointing  information  is  on  the  audio  track 


in  verbal  form  and  in  the  boresight  photographs.  Thus,  the  point  in  the 
ion  or  neutral  cloud  at  which  a  particular  spectral  line  shape  was  taken 
is  determined  by  a  study  of  the  Brush  recorder  charts  (with  real-time 
annotations  added  in  longhand),  the  audio  tracks  and  the  sequence  of 
boresight  photographs.  From  this  an  event  mark  is  added  to  the  l/2" 
magnetic  tape  and  an  entry  made  in  the  log  noting  time  and  pointing 
information.  As  a  result  of  the  inherent  inefficiency  of  this  process 
we  did  not  transcribe  other  Secede  III  events  to  1/2"  .magnetic  tape,  since 
other  tasks  of  the  contract  took  priority. 


Secede  III  Data  Anal 


Background  Information 


The  45  mm  aperture  Fabry-Perot  was  operated  from 
Ester  Dome  Observatory  near  Fairbanks,  Alaska  to  observe  the  Secede  III 
barium  releases  in  February  and  March  1969.  A  brief  survey  of  the 
purpose  of  the  releases  and  a  summary  of  cloud  morphology  will  serve 
to  introduce  the  problems  we  have  attempted  to  attack  in  the  data 
reduction. 

Historically  the  first  barium  releases  were  conducted  by  a 

12^ 

group  at  the  Max  Planck  Institute  for  Physics  and  Astrophysics  }  '  J 

for  providing  information  on  interplanetary  electric  and  magnetic  fields. 
Their  preliminary  tests  were  conducted  in  the  ionosphere  and  in  addition 
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to  evaluating  different  release  techniques,  gave  information  on  iono¬ 
spheric  electric  and  magnetic  fields  and  neutral  winds.  Subsequent 

4 

efforts  by  groups  in  the  United  States  have  been  directed  toward  using 
the  bariua  plasma  as  a  probe  to  test  ionospheric  properties  and  to 
attain  an  understanding  of  the  development  of  plasma  instabilities  in  the 
cloud.  For  these  studies  the  distribution  of  ions  in  space  roust  be 
determined.  Since  it  is  also  desirable  to  achieve  the  highest  possible 
release  efficiency  (ioniaable  bariua),  observations  of  the  evolution 
of  the  neutral  species  are  of  interest. 

The  bariua  releases  are  achieved  by  igniting  a  thermite-type 
reaction  in  cannisters  sent  aloft  in  carrier  rockets.  The  cannisters 
contain  Ba  and  CuO  in  a  molar  ratio  of  1.7:1*  The  reaction: 

1.7  Ba  +  CuO  -*  Cu  +  BaO  +  0.7  Ba  (+  96.3  kcal/mole)  ( 1  ; 

provides  neutral  atomic  bariua,  and  enough  energy  to  eject  a  spray 
of  superheated  liquid  reaction  products  from  the  cannister.  Neutral 
bariua  is  the  most  volatile  reaction  component  and  is  vaporized  from  the 
droplets.  The  first  ("snowplow"*)  phase  of  the  cloud  expansion  is 
extremely  rapid  (radial  expansion  velocities  on  the  order  of  1  km/sec) 
until  the  expanding  gas  has  swept  out  a  mass  of  ambient  significant 
compered  to  the  mass  of  released  gas.  This  takes  a  time  on  the  order  of 
0.5  sec.,  depending  on  the  size  of  the  release.  The  expansion  then 
slows  down  as  the  mass  of  the  swept-up  ambient  gas  becomes  large  compared 
to  the  mass  of  released  gas.  Is  can  be  expected  that  during  this  period 


*By  a  snowplow  expansion  we  imply  an  expanding  "solid"  sphere  which  drives 
the  ambient  atmosphere  outwurd. 
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the  ambient  atmosphere  (which  had  been  excluded  from  the  sphere  of 
expanding  reaction  products  during  the  first  phase)  will  begin  to  inter¬ 
mingle  with  the  expanding  barium  vapor.  The  second  phase  gradually 
degenerates  into  a  purely  diffusive  expansion  of  Ba  vapor  in  the  ambient 
atmosphere.  The  cooled  liquid  droplets  then  continue  outward  on  a 
ballistic  trajectory,  leaving  behind  a  sphere  of  neutral  bariun  vapor. 

The  expanding  particulate  sphere  has  no  further  effect  on  cloud 
morphology  after  the  atomic  bariun  has  been  vaporized.  The  release  must 
be  made  under  solar  illumination  for  an  ionized  cloud  to  be  produced, 
and  optical  observation  of  the  cloud  requires  a  reasonably  dark  sky 
background.  As  a  consequence  the  releases  are  done  during  late  twilight 
or  just  before  dawn.  In  the  absence  of  significant  cloud  drift  about 
twenty  minutes  of  observation  are  possible  in  either  case. 

Creation  of  the  major  portion  of  the  ionized  bariun  involves  a 
two-step  photoionization  process.  Figure  2  shows  partial  energy  level 
diagrams  for  neutral  and  ionized  barium  (Bal  and  fell) .  Ground  state 
neutrals  are  first  pumped  (by  solar  radiation  or  thermal  collisions)  into 

1  3 

the  D0  and  B  metastable  levels.  From  these  levels  they  are  then 
c.  df  x 

ionized  by  solar  radiation.  The  maximun  wavelength  for  ionization  of  the 

■j  G 

XD  metastable  is  X32&5A,  where  the  solar  flux  is  approximately  thirty 
3  0 

times  that  at  X2300A,  which  is  the  approximate  maximun  wavelength  which 
can  lead  to  ionization  of  ground  state  bariun.  A  prime  argument  for  the 
validity  of  the  two-step  photoionization  process  has  been  that  strontiun 
releases  provide  no  ion  cloud,  despite  the  fact  that  the  strontiun  term 

diagram  is  qualitatively  similar  to  bariun.  The  major  difference  is  that 

3  1 

the  'B  level  in  Sr  is  not  metastable,  and  the  D  metastable  does  not 


7 


ir 


appear  to  bo  pumped,  perhaps  because  of  its  proximity  in  energy  to  the 

1P  level/* 
o 

Creation  of  the  major  portion  of  the  ion  cloud  occurs  within  one 

k 

or  two  minutes  of  release.  Rosenberg  et  al.  give  a  time  constant  for 
ion  creation  at  170  km  of  close  to  10  seconds.  After  that,  however,  ions 
continue  to  be  created  frcm  the  neutral  cloud  at  a  slow  rate,  creating  a 
tenuous  "ion-bridge"  in  the  space  between  the  separating  ion  and  neutral 
clouds. 

Separation  of  the  main  clouds  is  caused  partly  by  neutral  winds 
carrying  the  neutral  cloud  away  frcm  the  ion  cloud  which  is  locked  to  the 
magnetic  field.  Also  any  electric  fields  which  are  present  will  result 
in  E  x  B  motion  of  the  ion  cloud  but  not  the  neutral  cloud.  The  neutral 
cloud  motion  has  been  used  for  mapping  winds,  and  the  ion  cloud  motion 
for  detecting  electric  fields,  especially  noticeable  in  the  auroral  zone. 
Mapping  of  density  in  the  clouds  by  standard  optical  techniques  (filter 
photography,  filter  photometers,  spectrcmeters)  is  greatly  complicated  by 
the  resonant  nature  of  the  optical  transitions  used  for  detecting  the 
cloud.  Trapping  of  this  resonance  radiation  destroys  the  simple  linear 
relation  between  intensity  and  line  of  sight  density  which  would  hold  in 
an  optically  thin  scatterer.  Consequently  estimates  of  total  inventory 
by  calibrated  photography  or  photometry  will  be  too  low. 

We  have  used  the  Fabj.y -Perot  to  measure  the  shapes  of  the  spectral 
lines.  These  measurements,  combined  with  calculations  of  h'jw  the  line 
shape  depends  on  line  of  sight  density,  give  a  much  better  picture  of 
actual  neutral  and  ionized  barium  density  distributions  in  the  clouds. 

The  problems  to  which  we  have  addressed  ourselves  are;  the  determination 
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of  Ba  I  ground  state  density  from  the  X5535  line,  important  for  cloud 
morphology  and  release  kinetics;  the  Ba  I  ("4) ?)  metastatic  density  from  the 
X5826  line,  bearing  on  ion  production;  and  the  Ba  II  ground  state  and 
metastable  state  densities  from  the  X4$)3^  and  X614-2  lines,  important  for 
characterizing  the  ion  cloud. 

The  calculations  of  resonance  radiation  transport  necessary  for 

interpreting  the  neutral  line  shapes  draw  heavily  on  the  theory  developed 

5 

by  Thomas  and  are  described  in  Appendix  A,  Briefly,  we  have  assumed  a 

Gaussian  radial  distribution  of  neutral  atoms,  expanding  only  by 

diffusion  with  a  diffusion  coefficient  of  0.05  km  /sec  (measured  by 

6 

AFCRL  on  the  Puerto  Rico  releases  at  188  lan).  Although  there  is 
considerable  hyperfine  structure  in  all  the  bariun  lines  due  to  isotope 
and  nuclear  spin  effects,  Table  I  shows  that  the  magnitude  of  the 
splitting  is  negligible  compared  to  doppler  temperature  broadening  except 
in  the  Ba  II  line.  Accordingly  for  the  X5535  radiation  transport 

calculations  we, have  assunecl  r.  single  frequency,  thermally  broadened 
line  shape  in  a  two  level  (ground  state  and  excited  state)  approximation. 
Further  approximations  are  discussed  in  Appendix  A. 

Line  shapes  were  obtained  on  all  six  releases  of  the  Alaska 
series,  but  data  from-  only  one  event  has  been  reduced  at  present.  This 
release,  event  "Gun",  is  of  interest  because  it  was  similar  in  size  and 
altitude  to  event  "Dogwood"  of  the  Puerto  Rico  series.  It  is  hoped 
that  comparison  of  the  two  events  will  be  of  interest  in  comparing 
differences  1 .  response  of  artificial  plasma  clouds  to  tne  arctic  and 
tropic  ionospneric  conditions. 
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TABLE  I 


Intensities  and  wavelength  shifts  for  byperfine  components  of  bariun  lines 
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Mass  5535  6142  4934 


Int 

Av(mK) 

Int 

£v(mK) 

Int 

Av(mK) 

138 

71.66 

0 

71.66 

0 

71.66 

0 

136 

7.81 

+  6 

7.81 

+  6 

7.81 

+  6 

134 

2.42 

+  9 

2.42 

+  9 

2.42 

+  9 

137 

1.87 

-1(5 

6.36 

-  5 

3.53 

-123 

3.74 

_  £ 

4.94 

+11 

3.53 

-  72 

5.61 

■(•16 

0.71 

+146 

3.53 

+197 

135 

1.10 

-10 

3.71 

-  5 

2.06 

-105 

2.20 

0 

2.89 

+  8 

2.0  6 

-  60 

3.30 

+23 

0.4l 

+136 

2.06 

+180 

NOTE  -  The  splitting  of  the  15535  line  would  not  be  completely  masked 
by  a  thermal  doppler  FWK4  of  33  mK  for  barium  at  900°K,  but  we 
would  not  expect  it  to  be  noticeaole  with  our  resolution  and 
signal-tc-noise,  as  was  the  case. 
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2.  Neutral  Cloud  -  X5535A  and  X  5826/1 

Figure  3  is  a  plot  of  the  optical  depth  effects  observed  in 
X >535  for  event  Gun.  As  in  studies  of  the  01  lines,  the  quantity  ve 
have  used  for  analysis  is  the  contrast  half -width  of  the  line  shapes,  i.e. 
the  full  width  at  an  intensity  halfway  between  the  maximun  and  minimun  in 
tue  profile.  In  determining  the  contrast  half -width  the  experimental 
points  have  been  i'ibteu  to  a  curve  using  the  first  half  of  tue  procedure 
detailed  in  Appendix  A.  The  fitted  curve  (before  deconvolution)  n us  then 
been  used  to  determine  the  contrast  finesse  (FSR  divided  by  contrast  haif- 
widtn).  The  family  of  curves  in  Figure  3  has  been  obtained  from  the 
resonance  radiation  transport  calculations  described  in  Appendix  in  Fach 
curve  represents  the  evolution  of  a  constant  inventory,  constant  tempera¬ 
ture  neutral  bariuu  cloud  in  time  with  diffusion  the  only  ’loss”  process. 
Instrumental  broadening  has  not  been  removed  from  the  experimental  hfuf- 


wJdths,  so  that  the  computer  generated  resonance  pro!  ij.es  nave  been 
convolved  with  the  appropriate  insurume..t  fvnetien  befeve  competing 
nuLf-width.  The  oomputacioixal  temperature  esod  for  the  family  of  curves 

in  Figure  3  is  5G0CK,  which  is  the  temperature  at  1'.  D  kn  given  by  too 

7 

USSA  (Supl. )  model  atmospnere  for  winter  and  ar.  exospheric  temperature 
ox  l‘.:tO°K  (obtained  from  tne  Jacchia  satellite  dra<.  model  mentioned 
earlier  with  the  parameters  fox*  Fairbanks  ac  Ob:  hO  Ur,  March  1 9, 

F10  y  =  1^0;  F10  Y  -  Kp  •-  1).  One  curve  was  also  run  for  T  -  ;%'CK 
to  show  how  the  inventory  would  change  if  a  different  cloud  temperature 
were  appropriate.  The  deduced  inventory  at  750°K  is  appro.cimuti../  ... 5 
times  the  inventor; y  a;  fGOlh.  The  deduced  inventory  for  a  900ca  cloud 
ut  midale  and  late  tines  is  constant  <;nu  equal  to  0. 3f>  of  the  total 
vaporizable  oariua  available  in  the  release. 
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If  the  neutral  barium  temperature  does  not  change  substantially. 

Figure  3  suggests  a  sharply  falling  neutral  inventory  in  the  early 

moments  of  the  cloud  history  when  ionization  of  the  neutrals  is  occurring. 

The  magnitude  of  this  fall  is  somewhat  deceiving.  If  the  initial  values 

are  taken  seriously  they  indicate  a  total  Ba  I  inventory  of  close  to 

7  x  10^  atoms,  compared  to  the  maxima  observed  yield  of  2C$,  or  1.3 
25 

x  10  atoms.  Therefore  either  the  transport  calculations  are  in  error 
or  some  of  the  observed  X5535  broadening  must  be  produced  by  the  tempera¬ 
ture  of  the  neutral  barium  being  substantially  in  excess  of  the  900°K 
ambient  atmospheric  temperature. 

We  have  considered  the  possible  distribution  of  the  excess  energy 
in  Reaction  (1  ),  _.e.  9 6  kcal/'mole.  Foppl  et  al.  require  60  kcal/mole 
to  bring  the  liquid  expelled  from  the  cannister  (reaction  products  plus 
excess  barium)  to  ~2000°K.  The  remaining  36  kcal/mole  is  then  available 
for  vaporizing  the  excess  (atomic)  barium, *  heating  the  vaporized  barium 
above  2000°K,  heating  the  entrained  and/or  penetrating  atmosphere,  and 
driving  the  initial  cloud  expansion.  If  we  assume  a  20$  efficiency  of 
vaporization  of  the  excess  barium  (the  maximum  according  to  measure¬ 
ment^  3  **  }  then  31  kcal/mole  remain  for  heating  and  expansion.  In  event 
Gum,  which  contained  W3  kg  of  reactants  with  composition  appropriate  to 
Eq.  (  1),  we  arrive  at  a  total  of  h.J  x  10^  kcal  (z  x  lO1^  ergs).  At  an 

<3  p 

ambient  pressure  of  3*5  x  10*°  dynes/cm  at  170  km  (USSA  Suppl., 

^exospheric  ~  1260°K,  winter  model)  ths  initial  "snowplow"  expansion  to  a 

13 

radius  of  1  km  requires  only  1. 5  x  10  ergs  to  displace  the  ambient 
atmosphere,  i.e.  less  than  10 $  of  the  available  energy. 


♦Latent  heat  of  vaporization  of  Ba  is  36  kcal/mole  at  1910°K. 
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The  energy  required  for  the  second  and  slower  expansion  phase 
(which  degenerates  into  the  diffusion  phase),  as  well  as  the  energy 
coasuned  in  heating  the  interdiffusing  ambient  atmospheric  species, 
depends  on  how  effective  the  hot  barium  vapor  is  in  continuing  to  drive 
back  the  ambient  atmosphere.  Models  of  the  early  expansion  period  are 
under  investigation  by  AFCRL  (Reference  8,  p.  5)  and  W.  B.  Murcray  of 
the  University  of  Alaska  (Reference  8,  p.  97). 

We  have  been  discouraged  fraa  pursuing  further  cloud  calculations 
on  the  basis  of  simplified  models  of  pure  adiabatic  snowplow  expansion 
or  complete  raix.\/ig  for  the  following  reasons.  Adiabatic  expansion  with 
no  mixing  could  enlarge  the  cloud  only  to  a  radius  of  2. 5  bn  with  the 
energy  available  (2  x  XO1^  ergs),  after  which  the  expansion  would  be  by 
diffusion  and  would  increase  the  cloud  radius  relatively  slowly 
(~200  m/sec).  However,  the  cloud  is  observed  to  grow  to  a  5  km  radius 
in  the  first  5  seconds  before  diffusion  takes  over.  Alternatively,  a 
complete  mixing  model  would  require  complete  sharing  of  the  excess  energy 
with  all  ambient  species  in  the  initial  5  km  sphere.  This  would  result 
in  a  temperature  rise  of  only  ~50°K.  But  there  is  certain  to  be  seme  pre¬ 
heating  of  the  surrounding  atmosphere  by  the  outgoing  liquid  droplets, 
and  seme  interior  volune  of  the  sphere  is  certain  tc-  be  evacuated  by  the 
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initial  (fast)  expansion.  Thus  neither  simplified  model  appears  to  he 
refined  enough  to  even  estimate  what  the  barium  temperature  should  be 
after  the  cloud  has  expanded  to  5  km. 

Common  sense  would  suggest  that  at  least  the  interior  of  the 
cloud  would  remain  at  an  elevated  temperature  for  some  time  after  initial 
expansion.  If  this  is  the  case  then  the  narrowing  of  the  experimental 
15535  lines  with  time  would  involve  both  a  depletion  of  neutral  inventory 
(through  ionization  and  oxidation)  and  a  cooling  of  the  cloud  (by 
conduction  and  intermixing).  The  time  constant  for  conductive  cooling  of 
a  1  km  sphere  of  Ba  and  Ba+  at  2500°K  grading  to  a  900°K.  ambient  tempera¬ 
ture  at  5  km  is  initially  between  10  and  100  seconds  depending  on  density 
inside  the  sphere  (lO^cnT^  -  10®  The  time  constant  for  diffusive 

intermixing  should  be  of  the  seme  order  of  magnitude,  based  on  the  observa¬ 
tion  that  the  Ba  cloud  expands  radially  at  roughly  200  ra/sec  when  it  is 
5  km  in  radius.  Whether  a  refined  cloud  expansion  model  could  allow  as 
high  an  initial  temperature  as  we  seem  to  observe  must  await  the  results 
of  the  calculations  mentioned  above. 

Figure  4  1  is  an  attempt  to  assess  the  competing  effects  of 
resonance  absorption  and  atom  temperature  on  the  line  profiles.  In  it  an 
experimental  line  shape  obtained  45  seconds  after  release  has  been 
compared  with  two  computer  generated  profiles  having  the  same  half-widths 
but  different  combinations  of  temperature  and  neutral  density.  The 
experimental  curve  has  been  folded  about  its  midpoint  to  make  the  distribu¬ 
tion  of  experimental  points  more  continuous.  The  computed  profiles  are 
appropriate  to  a  cloud  of  5  km  initial  radius  (measured  value  appropriate 
to  R  +  5  sec.)  having  undergone  45  seconds  of  diffusive  spreading.  It  is 
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seen  that  a  very  good  fit  is  provided  by  a  curve  with  a  ground  state 
inventory  of  2.8  x  10  (which  is  near  the  value  indicated  in  Figure  3 
for  late  times)  and  a  temperature  of  2000 °K,  substantially  above  the 
expected  ambient  temperature  of  900°K  for  this  altitude.  It  thus 
appears  that  by  very  early  times  the  bulk  of  neutral  depletion,  i.e. 
ionization  end  oxidation,  may  have  taken  place,  and  that  subsequent 
ionization  produces  a  negligible  drain  on  the  neutral  inventory.  The 
rapid  decrease  of  line  width  indicated  in  Figure  3  seems  to  be  caused 
almost  entirely  by  cooling  processes.  The  fact  that,  after  5  minutes,  the 
data  follow  a  constant  temperature,  constant  inventory  model  curve 
suggests  that,  if  continued  conversion  of  Ba  to  BaO  is  taking  place,  it 
must  be  accompanied  by  a  back  reaction  which  replaces  free  bariun  as 
fast  as  it  is  lost.  It  should  be  noted  that  the  effect  of  using 
different  instrumental  defect  functions  in  generating  the  predicted 
curves  in  Figure  4  has  been  investigated  and  found  to  be  insignificant. 

The  physical  parameter  which  is  most  questionable  in  these  cal¬ 
culations  is  the  temperature.  Figures  3  and  4  have  shown  that  the 
inferred  neutral  inventory  is  very  sensitive  to  the  choice  of  temperature 
of  the  radiating  atoms.  It  can  be  safely  assuned  that  at  late  times 
the  cloud  has  reached  the  temperature  of  the  ambient  atmosphere,  so  that 
the  constancy  of  the  neutral  inventory  at  late  times  is  Justified.  One 
way  of  determining  the  temperature  independently  would  be  from  the 
rotational  structure  of  the  BaO  bands  observed  in  the  release.  Since  in 
another  experiment  a  spectrometer  was  observing  these  bands  (Reference  8, 
p.  86)  during  the  Gun  event,  it  may  be  possible  to  obtain  the  desired 
information  in  the  near  future. 
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Such  temperature  measurements  from  BaO  hands  have  been  reported 
by  Harang  and  Stoff regent  from  Ba  releases  by  the  Max-Flanck-Institute. 
They  recorded  BaO  temperatures  at  R+l  minute  appropriate  to  undisturbed 
atmospheric  temperatures.  They  do  not  specify  the  size  of  the  releases, 
but  typically  the  Max-Planck  releases  are  smaller  than  Gun  by  a  factor 
of  50.  Such  a  smull  release  would  tend  to  equilibrate  much  faster,  so 
that  the  Harang  and  Stoffregen  measurements  do  not  necessarily  apply  to 
early  times  of  the  large  Alaska  releases. 

Without  knowing  the  distribution  of  BaO  within  the  cloud,  it  would 
be  hazardous  to  apply  BaO  temperature  measurements  to  the  question  of 
whether  the  cloud  temperature  is  high  enough  to  thermally  populate  the 
metastahle  levels.  If  the  BaO  is  generated  by  the  ambient  atmosphere  as 
it  penetrates  and  cools  the  outer  reaches  of  the  cloud,  then  the  BaO  tem¬ 
peratures  would  be  much  lower  than  the  interior  temperature  of  the  cloud 
where  ion  generation  is  taking  place.  Thus  a  low  BaO  temperature  would 
not  necessarily  contradict  a  high  interior  temperature,  whereas  a  high 
BaO  temperature  would  support  it.  The  BaO  measurements  can  therefore 
only  be  treated  as  a  lower  limit  to  the  interior  cloud  temperature. 

Another  method  of  determining  the  cloud  temperature  would  be  to 
determine  the  width  of  the  X5826  metastable  line  at  very  low  intensity 
so  that  optical  depth  effects  are  negligible.  This  method  has  the 
advantage  that  it  gives  the  temperature  at  the  location  in  the  cloud 
where  the  metastable  level  is  being  populated.  Unfortunately  the  profiles 
we  obtained  of  the  X5826  line  are  entirely  too  noisy  for  this  purpose. 

Instead  of  determining  the  exospheric  temperature  ss  we  have  done 
for  Gum  by  using  measured  decimeter  flux  and  magnetic  index  values  in 


conjunction  with  Jacchia's  satellite  drag  model,  it  would  be  more 
accurate  to  use  the  interferometer  to  measure  the  exospheric  temperature 
from  the  01  X6300  line  sometime  during  the  course  of  the  cloud  evolution. 
In  either  case  the  exospheric  temperature  must  be  extrapolated  to  release 
altitude  through  a  model  atmosphere,  but  the  measured  red  line  tempera¬ 
ture  would  remove  at  least  one  uncertainty.  In  practice,  accumulating  a 
reasonable  red  line  profile  would  probably  take  sufficiently  long  that 
it  would  be  preferable  to  accumulate  the  red  line  after  the  barium  cloud 
has  disappeared  and  then  extrapolate  backward  in  time  using  the  Jacchia 
model. 

One  ;  comparison  that  has  proven  interesting  in  all  the  barium  lines 
is  between  the  line  intensity  at  the  peak  of  the  Fabry-Perot  profile,  the 
integrated  FP  profile  intensity  and  the  integrated  line  intensity  which 
comes  through  the  monitor.  In  Figure  5  the  X5535  data  points  show  the 
measured  time -.variation  of  these '.quantities.  AIbo  shown  are  the  varia¬ 
tions  predicted  by  the  computer  solution  for  a  900°K  diffusing  cloud  of 
constant  inventory  2  x  10  J  atoms.  Two  things  are  especially  striking 
about  this  plot.  The  first  is  the  rapid  initial  decline  of  the  monitor 
intensity  which  is  not  matched  by  the  FP  peak  intensity.  Our  first 
impression  was  that  this  resulted  from  contamination  by  BaO  bands. 

Brown  (Ref.  7,  p.  86)  has  detected  an  interfering  band  in  a  spectrum 
> 

taken  during  the  first  40  seconds  of  Gun,  but  published  no  intensities. 
When  the  FP  profile  is  integrated  over  the  line,  however,  the  square 
points  in  Figure  5  result,  indicating  that  no  major  contamination  by 
BaO  is  present.  It  is  our  present  feeling  that  the  initial  rapid  fall 
of  the  monitor  and  integrated  FP  X5535  intensities  is  linked  to  the  high 
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initial  temperature  during  the  creation  of  the  hulk  of  the  ions. 
Although  analyses  similar  to  Figure  4  will  he  applied  to  all  the  early 
profiles  in  an  attempt  to  separate  the  cooling  effect  from  the  ionization 
effect,  it  is  not  expected  that  these  analyses  will  he  too  successful. 

The  reason  is  apparent  frcm  looking  at  the  scatter  of  points  near  the  peak 
of  the  profxxe  in  Figure  4„ >  Much  depends  on  the  placement  of  the  peak 
intensity,  and  the  scatter  is  too  large  to  do  this  with  any  great  certain¬ 
ty.  In  future  releases  we  will  attempt  to  obtain  a  better  early  time 
history  of  the  X5535  line  profiles. 

Even  with  the  data  in  preliminary  form,  however,  a  question  is 
raised  with  respect  to  the  two-step  photoionization.  Since  the  Bal 
inventory  remains  stable  for  a  long  period  of  time  even  under  solar 
illmination  as  intense  as  that  flux  responsible  for  the  initial  ioniza¬ 
tion  and  since  the  main  ionizf'.oion  is  produced  initially  when  the  plasma 

is  very  hot,  i;  appears  that  cue  hi  ■c;  tore  of  the  barium  cloud  is 

2.  3 

more  efficient  in  populating  the  f  C  metastable  levels  than  solar 

pumping.  In  face  if  thermal  equilibrium  were  established  at  2500*%  a 

3 

reasonable  temperature  in  very  early  time,  the  D  levels  which  lie  1.13 
eV  above  the  ground  state  and  together  have  a  statistical  weight  of  15 
would  contain  T .  37*  of  the  total  neutral  inventory,  and  the  "4)  level  with 
statistical  weight  5  which  lies  1  44  eV  above  the  ground  state  would 
contain  0.6^  of  the  neutral  inventory.  To  our  knowledge  the  themal 
population  of  metastable  levels  has  not  been  considered  previously  as 
the  major  source  of  metastables  which  are  subsequently  pbotoionized. 

The  second  striking  feature  of  Figure  5  is  the  continued  and 
rapid  decline  of  intensity  of  the  monitor  and  FP  peak  at  times  when  the 
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cloud  model  predicts  much  less  of  a  decline  for  both 
of  those  quantities.  If  the  cloud  maintains  a  constant  neutral 
bariua  inventory  throughout  its  diffusive  expansion  (as  the  recorded 
profiles  suggest),  then  even  at  R  +  18  minutes  (0537)  the  cloud  has  an 
optical  thickness  of  2.5  and  is  scattering  all  the  sunlight  incident  upon 
it.  The  model  then  predicts  a  very  modest  decline  in  peak  intensity.  If 
our  interpretation  of  the  profiles  is  correct,  then  the  falling  intensity 
must  result  from  a  decline  in  the  exciting  solar  flux  caused  by  atmo¬ 
spheric  scattering.  Since  the  last  data  point  at  R  +  l8  minutes  was 
taken  just  before  the  neutral  cloud  vanished  (when  the  earth’s  shadow 
height  reached  170  km)  this  interpretation  seems  reasonable  for  very 
late  times.  On  the  other  hand,  whether  one  would  expect  to  find  a 
decreasing  solar  flux  irradiating  the  cloud  at  early  times  must  await 
more  accurate  determinations  of  the  cloud  height  frcm  triangulation,  as 
well  as  numerical  calculations  of  atmospheric  scattering.  Results  of  a 
detailed  calculation  of  atmospheric  attenuation  have  been  published  by 
Lloyd'*'®  and  will  be  applied  when  cloud  position  data  become  available. 

A  crude  graphical  calculation  (ignoring  refraction  effects)  shows  that 
at  release  time  the  solar  flux  illuninatihg  the  cloud  passes  only  40  km 
above  the  surface  of  the  earth  and  tnis  height  decreases  by  ~2  km  per 
minute. 

The  observation  of  steady  neutral  Ba  inventory  despite  falling 
line  intensity  from  the  cloud  exhibits  another  shortcoming  of  photo¬ 
metric  observations,  i.e.  that  in  de /eloping  a  time  history  of  the  cloud 
the  time  dependence  of  the  illuminating  intensity  must  be  known,  whereas 
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the  density  determinations  from  line  shape  observations  are  insensitive 
to  the  intensity  of  the  flux. 

Figure  6  shows  the  time  variation  of  the  FP  and  monitor 
intensities  of  the  X5826A  line  which  has  the  metastable  state  as  a 
lower  level.  The  abrupt  change  in  slope  coincides  in  time  with  the 
break  observed  for  the  X5535  line.  The  appearance  of  the  slope  transi¬ 
tion  in  both  wavelengths  supports  our  belief  that  there  are  two 
separate  processes  involved  in  populating  the  metastable  levels.  In  our 
model  the  break  in  the  intensity  curves  would  coincide  in  time  with 
a  change  in  the  process  responsible  for  pumping  the  metastable  levels. 

We  would  expect  a  corresponding  break  in  the  intensity  curves  for  lines 
terminating  on  the  D  levels,  but  occurring  later  in  time  than  for  the 
X5826  line  since  the  Boltzman  factor  is  less  for  the  lower  lying  "D 
levels  and  they  should  remain  in  thermal  equilibrium  longer  with  the 
ground  state.  Since  the  line  peak  intensity  tracks  the  integrated 
(monitor)  intensity  so  well,  it  appears  that  at  no  time  is  the  bariun 

cloud  optically  thick  for  X5826.  An  upper  limit  may  thus  be  placed  on 

1  1 
the  D  density  and  translated  into  percentage  of  neutrals  in  the  D 

state  for  given  ground  state  densities.  At  R  +  45  seconds  this  limit 

is  either  <3$  or  <60 *[*,  using  the  ground  itste  densities  deduced  from 

Figure  18  at  cloud  temperatures  of  900°K  or  200O°K,  respectively.  The 

difference  of  slopes  at  late  times  and  low  signal  levels  is  assuned  to  be 

the  result  of  contamination. 

Pictures  taken  with  the  boresight  camera  snowed  that,  while  we 
were  generally  not  pointed  at  the  exact  center  of  the  cloud,  we  were  less 
than  one-quarter  of  a  Gaussian  radius  away  frcm  the  center.  Calculations 
of  the  distortion  of  the  X5535  line  profile  generated  by  this  sighting 
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error  reveal  that  within  one-quarter  of  the  Gaussian  radius  the  half¬ 
width  should  change  less  than  3^-  Doubtless  some  of  the  fluctuation  of 
the  experimental  points  in  Figure  3  is  due  to  sighting  closer  or  farther 
away  from  the  cloud  center.  Since  boresight  pictures  were  not  taken 
continuously  enough  to  provide  a  complete,-  record  (l  picture  every  3-4 
scans  at  early  times)  we  have  ignored  the  small  sighting  errors  and 
assumed  in  all  calculations  that  we  were  aimed  at  the  cloud  center. 

3.  Ion  Cloud  -  X4934A  and  X6l42A 

In  addition  to  the  Bal  lines  X5826  and  X5535>  measurements  were 
carried  out  on  the  X4934  and  X6l42  lines  of  Ball.  The  X6l42A  line 
terminates  on  the  nietastable  state  of  Ba+  and  we  had  planned  to 

obtain  metastable  densities  by  observing  optical  depth  effects. 

Unfortunately  the  FP  spacer  was  adjusted  for  the  Gun  release  so  that  in 
scanning,  less  than  half  a  X6l42  profile  was  recorded,  making  interpreta¬ 
tion  impossible.  The  X6l42  monitor  intensity  was  comparable  to  X4934, 
so  that  density  estiuates  should  be  possible  on  other  releases  of  the 
series  where  the  X6l42  line  placement  in  the  FP  scan  was  more  favorable. 

As  shown  in  Table  I  (p.  ?>)  the  X4934  ion  line  is  the  only  one  of 
the  four  lines  observed  where  the  hyperfine  splitting  is  large  enough  to 
noticeably  affect  the  line  shape.  Using  the  values  of  the  splittings  and 
relative  intensities  given  in  Table  X  (obtained  from  Kivel  et  al. a 
line  shape  was  generated  for  an  optically  thin  X4934  source  at  800°K*  and 
convolved  with  an  instrumental  function  appropriate  to  late  time  observations 


*The  choice  of  800°K  for  the  X4934  computational  ion  temperature  was  an 
estimate  made  during  preliminary  data  reduction.  The  heuristic  nature  of 
Figure  7  does  not  justify  the  effort  of  recomputing  the  profile  and  re¬ 
drawing  the  figur."  with  a  more  accurate  temperature  of  900°K. 
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of  the  ion  cloud.  This  profile  is  shown  in  Figure  7  along  .» it'i  a 
>rofile  observed  at  R  +  26  minutes  in  Gun.  A  complete  profile  extends 
from  the  values  52  to  173  on  the  arbitrary  frequency  scale;  partial 
profiles  from  the  interference  order  before  and  after  the  complete  profile 
are  also  shown.  The  eleven  hyperfine  components  have  been  combined  into 
J  groups  of  closely  spaced  lines  v  .th  centers  of  gravity  indicated  by  the 
arrows  labeled  with  each  group's  relative  intensity.  The  thermal  doppler 
broadened  jrofile  of  the  isolated  component  on  the  left  of  the  main 
component  is  decoupled  from  radiation  in  the  main  peak  and  will  henceforth 
be  referred  to  as  the  "hyperfine  line".  The  experimental  shape  has  been 
normalized  to  the  peak  of  the  hyperfine  line  and  the  minimum  to  the  left 
of  the  hyperfine  line  of  the  theoretical  profile.  It  is  obvious  that 
even  at  such  1’te  times  in  the  cloud  history  there  is  substantial  absorp¬ 
tion  of  the  main  component,  and  that  such  absorption  will  lead  to 
erroneous  estimates  of  line  of  sight  ion  densities  from  photometry  aid 
photometric  photography.  For  an  optically  thin  medium  the  ratio  of  cne 
intensify  of  t^e  main  peak  to  the  hyperfine  peak  should  be  ih,7.  -it  late 


times  this  ratio  approached  10  in  the  ion  cloud,  indicati.ig  an  optical 

depth  of  about  0.9  for  the  center  of  the  nain  component,  or  a  i:.ne  of 

Xl.  p 

s  :l  _;ht  _ens:lty  of  about  1,5  x  10  ground  state  ions  per  cm  (see  Appendix 
A  for  re  .at ‘on  between  main  no  hyperfine  ratio  ar.d  optical  depth'. 

Figure  8  shows  the  time  variation  of  the  monitor  (integrated 
total  intensity),  the  main  peak,  and  the  hyperfine  peak  intensity  of  the 
1*034  line.  It  appears  that  the  hyperfine  intensity  rather  quickly 
attains  a  steady  (exponential)  decay  rate  which  is  only  reaches  at  a  much 
,‘.l> tor  time  ‘  y  the  main  peak  intensity  hince  our  earlier  resonance 


Reproduced  from 
best  available  copy. 


22 


radiation  calculations  have  shown  that  an  optically  thick  scatterer 
should  give  a  reasonably  constant  intensity,  it  is  our  impression  that 
the  main  component  decay  must  be  largely  due  to  the  falling  solar  flux 
intensity. 

The  large  scatter  of  data  points  around  the  fitted  curves  in 
Figure  8  is  caused  by  slight  pointing  changes  with  respect  to  the  fine 
striations  in  the  ion  cloud.  At  about  R  +  6  min.  for  Gun  the  rf  returns 
underwent  a  change  in  nature  (Ref.  8,  p«  175)  from  what  would  be 
expected  of  a  single  scattering  object  to  a  multiple  object.  This  time 
coincides  reasonably  with  the  visual  change  of  the  ion  cloud  from  a 
smooth  to  a  striated  shape;  looking  eventually  not  unlike  a  rayed 
auroral  arc,  but  lacking  the  rapic  internal  motion  characteristic  of  a 
ra.yec  ere,-  Sarly  in  time  the  ion  cloud  may  be  characterised  rather 
easily  by  a  diffusion  model  with  different  diffusion  coefficients 
parallel  and  perpendicular  to  the  magnetic  field  and  ion  density  ue- 

creesing  mor.ctonicsJ.ly.  Prior  to  stristion  the  ion  density  appears  to 

12 

begin  to  increase  again  until  the  plasma  instability  ta.;e.c  over 
completely  and  partitions  the  ior.  cloud  into  a  multitude  of  filamentary 
bundles.  For  Gum  trie  F?  was  pointed  at  the  region  of  peak  ion  light 
intensity  from  roughly  to  seconds  before  striation  (as  determined  by  the 
appearance  of  multiple  returns  from  the  SRI  short  pulse  radar 
(Rif.  7*  p.  IT.-  ))  tuatil.  some  minutes  later.  The  individual  scans  of 

will  be  investigated  very  thoroughly  when  complete  reduction  of  the 
bariuc  data  is  undertaken  to  obtain  ion  colunn  densities  nr -.or  to  and 
Just  after  stritlon.  The  chief  difficulty  with  analysis  of  the  Fabry - 
Perot  data  is  that  we  have  no  idea  of  the  scale  of  the  structure  in  the 
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line  of  sight,  so  that  colunu  densities  may  not  be  easily  converted  to 
volume  densities,  the  quantity  of  interest  for  comparison  with  rf 
measurements.  This  problem  is  particularly  bothersome  after  striation, 
when  the  scale  of  the  structure  may  be  even  smaller  than  the  acceptance 
cone  ct  the  F?.  Integration  over  several  filaments  which  have  voids 
between  would  give  us  a  lower  estimate  of  ion  densities  than  would 
actually  be  present. 

The  problem  of  estimating  peak  densities  and  structure  scale  in 

the  striated  icr.  cloud  is  now  under  attack  using  the  hyperfire  line.  A 

new  Fabry- ’>erot  is  under  construction  for  use  on  the  Secede  II 

barium  releases  which  will  be  matched  to  a  k-power  fore  telescope  to  give 

k  minutes  of  arc  ra  rad.)  angular  resolution.  This  is  equivalent  to 

100  meter  resolution  at  a  slant  range  cf  200  km.  With  the  chamber 

pressure  fixed  at  a  point  in  the  scan  corresponding  to  the  peak  in  the 

hyperfine  line,  a  spatial  scan  will  la  executed  across  the  strirtions. 

This  will  tell  us  the  ion  cloud  fine  scale  striation  structure  via  a  .line 

,  th 

with ai  absorption  coefficient  per  atom  which  is  1/15  that  of  the  main 
line . 

If  we  arc  willing  to  consider  e  20f>  deviation  from  the  linear 
relation  between  th-  h.f.  line  intensity  and  colrnn  density  as  referring 
to  a  r.odiir.  -hict  is  "optically  thin",  then  it  should  be  pcs-'cle  to 
quantitatively  characterize  the  beriun  lor.  cloud  at  points  Jr.  its  time 
history  "her.  the  ms-’r  to  hyperfir.e  ratio  is  as  low  as  According  to 

the  fitted  our-'is  in  Fig.  7  this  occur:-:  at  R  <-  ”  minutes,  ..  ...>rt  a  after 
striation.  However,  there  is  a  single  pc-.’nt,  which,  may  re  -r. sent  a  look 
at  the  nidi]  -)  of  a  striation,  at  R  ■?  1C  minutes  where  the  ratio  is  less 
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than  2.5*  Thus  we  may  have  to  take  care  in  interpreting  the  spatial 
scans.  With  the  instrument  calibrated  for  absolute  intensity  and  inter¬ 
spersing  spatial  scans  with  spectral  scans  to  tell  us  what  the  main/ 
hyperfine  ratio  is  for  a  giver  hyperfine  intensity,  we  should  be  able  to 
make  approximate  corrections  for  optical  depth  in  the  hyperfine  line. 

For  a  release  the  size  of  Gum,  the  hyperfine  line  should  rot  exceed 
optical  depth  unity  except  for  very  early  times  in  the  formation  and 
diffusive  expansion  period. 

The  farthest  we  can  presently  take  the  analysis  of  the  Alaska 
data  is  to  obtain  an  estimate  of  peak  column  densities  of  ground  state 
ions.  This  comes  most  easily  from  the  use  of  Figure  U  in  Appendix  A 
relating  to  the  interpretation  of  the  main/hyperfine  ratio.  Tae  meta¬ 
stable  population  must  be  considered  also;  Duell  and  Sears, ^  from  a 
consideration  of  transition  probabilities  and  solar  flux,  estimate  that 
6C '/»  of  the  ion  population  is  in  the  ^/g  metastable  levels  (c.f. 

Figure  2).  Furtheimore  the  ion  clouds  are  aligned  with  the  magnetic 
field  and  hence  are  aspect  sensitive.  When  photographic  determinations 
of  the  size  of  individual  striations  become  available,  we  will  attempt 
to  generate  ion  densities  to  compare  with  rf  measured  electron  densities. 
This  must  be  done  with  extreme  care  because  in  addition  to  observing  the 
striations  from  an  aspect  other  than  normal  to  the  field  lines,  we  could 
be  looking  through  several  striations,  or  at  some  striations  which  are 
partially  shielded  from  the  sun  by  others,  with  a  resultant  distortion 
of  the  uniform  solar  flux.  If  the  ion  cloud  were  a  picket  fence,  analysis 
would  be  straightforward.  Photographs  looking  up  the  magnetic  field 
lines  show  that  it  has  more  the  appearance  of  an  auroral  corona.  For 
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these  reasons,  obtaining  ion  densities  from  optical  depth  measurements 
is  very  hazardous,  and  has  not  yet  been  attempted  in  data  reduction  for 
the  Alaska  releases. 

Alternatively  it  is  possible  to  go  in  the  other  direction  and  to 

use  the  rf  ng  measurements  to  derive  an  effective  structure  size  which 

would  apply  if  we  were  looking  at  only  a  single  uniformly-illuminated 

striation.  For  instance  at  R  +  (10  min.,  20  sec.)  the  main/hyperfine 

12 

ratio  was  2.5>  implying  a  ground  state  column  density  of  1.3  x  10 

ions/cm2.  This  gives  a  total  (ground  state  plus  metastable)  ion  column 

12  -2 

density  of  3.2  x  10  cm  .  Considering  that  our  viewing  aspect  to  the 
field  lines  was  about  k0°  (estimated  at  present;  more  accuracy  will  come 

from  triangulated  cloud  positions),  the  ion  column  density  normal  to  the 

12  -P  7 

field  lines  is  2  x  10  cm  .  At  this  time  the  rf  n  values  are 

3.3  ,  rf  .-3  H3  and  W>,  ,B  ,  rf  J  <«,  _  ateP 
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sounder)  and  1.1  x  10  cm~J  (SRI,  Kodiak  step  sounder).  Thus  the  single 
structure  size  that  this  type  of  analysis  would  indicate  at  R  +10 

12  6  e 

min.  is  ^  x  10  /k  x  10  =  5  x  Kr  cm  or  5  ton-  This  seems  extremely 
large  for  a  single  striation  and  must  represent  the  sun  of  several 
striations. 

The  final  type  of  information  which  can  be  obtained  from  FP 

profiles  is  drift  velocities  of  the  various  clouds.  The  only  occasion 

we  have  had  to  examine  our  Alaska  data  for  such  effects  was  when  it  was 

14 

mentioned  that  the  event  "Elm"  ion  cloud  appeared  to  be  moving  west¬ 
ward  at  supersonic  velocities.  From  a  rough  look  at  the  real-time  chart- 
records  we  confirmed  the  magnitude  and  time  dependence  of  velocity  for 
the  Elm  ion  .cloud  as  deduced  from  the  University  of  Alaska  TV  tracking 


system.  Velocity  information  is  inherent  in  all  our  records  of  both 
neutral  and  ion  cloud,  with  a  velocity  resolution  of  ~  ±30  meters/second 
for  the  better  line  profiles.  These  drift  velocity  determinations  give 
information  on  electric  fields  and  neutral  winds.  But  since  cloud  drift 
can  so  easily  be  obtained  from  the  TV  tracking  or  all  sky  camera  systems, 
and  since  we  have  many  other  problems  which  only  the  FP  can  solve,  we 
probably  will  not  investigate  the  velocity  data  further. 
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III.  Secede  II  Apparatus  Construction,  Field  Observations,  and  Initial 
Data  Reduction 
A.  Apparatus 

As  noted  in  the  previous  section,  a  larger  (100  mm  aperture) 
Fabry  Perot  interferometer  in  a  high-stability  mounting  cage  was  to  be 
fabricated  for  use  in  conjunction  with  a  16"  (kX  magnification)  telescope 
to  provide  high  spectral  resolution  at  higher  optical  speed  and  better 
spatial  resolution  (l/2  mrad,  or  100  meters  at  200  km  slant  range).  Also, 
electronic  data  processing  to  provide  real-time  transcription  of  the  data 
onto  l/2"  magnetic  tape  in  a  format  which  could  be  readily  reduced  with 
computer  assistance  was  a  requirement  cf  the  new  apparatus,  The  Fabry- 
Perot  optical  interferometric  system  and  the  electronic  control  system  for 
the  instrument  were  developed  by  our  group  at  the  University  of  Pittsburgh; 
the  16"  telescope  pointing  system  for  the  Fabry-Perot,  bore-sight  TV  and 
videotape  recorder,  and  the  electronic  data  transcription  systems  were 
developed  by  General  Electric  Space  Science  Laboratories,  Valley  Forge,  Pa. 

In  order  to  maintain  parallel  alignment  (to  ~  l/lOO  wavelength)  of 
the  Fabry-Perot  plates  during  the  inevitable  accelerations  and  vibrations 
associated  with  operation  on  the  telescope  training  mount,  a  special  invar 
mounting  cage  was  constructed  which  employed  controlled  deformation  of 
very  rigid  structures  to  achieve  and  maintain  plate  parallelism.  This 
design  was  modified  somewhat  from  the  original  version  in  order  to  improve 
the  smoothness  of  the  parallelism  adjustment.  In  addition,  a  wavelength 
sensitive  beam- splitting  mirror  and  dual  detection  head  was  attached  to  the 
output  stage  of  the  interferometer  to  permit  simultaneous  scanning  of  two 
spectral  lines  (as  opposed  to  the  sequential  line  scanning  system  of  the 
Alaska  instrument)  without  loss  of  light  in  either  channel.  In  this  instru¬ 
ment,  a  separate,  boresighted  monitor  filter  photometer  was  added  to  provide 
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total  intensity  records  of  the  various  lines  (X4934,  X5575>  X6lll,  X6l42) 
observed  by  the  interferometer. 

An  electronic  control  unit  making  extensive  use  of  integrated 
circuits  was  constructed  to  permit  flexible  operation  of  the  interferometer 
in  the  wavelength  scan  and  wavelength-lock  (spatial  scan)  modes,  e.g, 
variable  rate  of  scan,  dwell-time  per  channel,  number  of  channels  per  scan . 
In  addition,  remote  control  of  all  changes  at  the  interferometer  (such 
as  filter  selection;  X5535  & X4934  in  one  channel;  X6328  [laser  calibrator], 
\6l42  and  X6lll  in  the  other)  was  provided.  The  output  pulses  of  the 
three  detecting  photomultipliers  were  fed  to  specially  designed  preamplifier, 
discriminator,  wave  shaping  and  cable  driving  circuits  to  provide  standard 
pulses  for  the  data  processing  circuits  built  by  G.E. 

B.  Field  Observations 

The  final  integration  of  the  Fabry-Perot  with  the  telescope 
augmented  pointing  system  took  place  at  the  field  site  S-3a  (30°  6'N, 

85 0  14-2'W)  at  Wewahitchka,  Fla.  The  optical  interferometer  systems 
functioned  properly,  achieving  a  detection  sensitivity  very  close  to  that 
predicted  for  the  Secede  II  releases  on  the  basis  of  our  Secede  III 
Alaskan  observations.  Most  of  the  difficulties  in  the  field  occurred  with 
the  spatial  scan  computer  which  controlled  the  mount  pointing,  with  the 
boresight  TV/videotape  system,  and  the  digital  magnetic  tape  data  recording 
system. 

These  deficiencies  were  corrected  sufficiently  well  by  the  time  of 
first  release  that  good  data  was  obtained  on  all  twilight-nighttime  releases 
(with  the  exception  of  Nutmeg,  where  only  the  real-time  Brush  recorder 
records  were  obtained,  owing  to  failure  of  the  digital  magnetic  tape  data 
system;  the  videotape  recorder  failed  on  both  Nutmeg  and  ?.1um).  The 
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rationale  behind  the  Secede  II  field  observations,  together  with  a  preliminary- 
analysis  of  some  of  the  data,  is  given  in  the  next  section. 

C.  Initial  Secede  II  Data  Reduction 

1.  Introduction 

Observation  with  Fabry-Perot  (FP)  interferometers  on 
barium  releases  in  Alaska1^  and  Puerto  Rico-1 3  have  shown  that  the  \4934  line 
profile  is  severely  distorted  by  optical  depth  effects.  These  effects  give 
rise  to  a  relationship  between  source  brightness  and  line-of-sight  ion 
density  which  is  unknown  because  of  non-unifom  source  excitation.  This 
means  that  photographs  taken  in  the  light  of  such  an  optically  thick  transi¬ 
tion  do  not  provide  reliable  quantitative  information  on  absolute  or 
relative  ion  column  densities. 

Line  profiles  from  the  Alaskan  series^  suggested  that  a  subsidiary 
component  of  the  \4934  line  which  is  well  removed  from  the  main  peak 
becomes  optically  thin  quite  soon  after  release  (1-3  minutes).  This  sub¬ 
sidiary  component,  hereafter  referred  to  as  the  hyperfine  (hf)  line,  was 
expected  to  be  uniformly  excited  in  the  cloud  and  provide  reliable  quanti¬ 
tative  information  on  ion  column  densities.  2fcus  an  instrument  (described 
in  the  preceding  section)  was  constructed  which  was  capable  of  isolating 
this  component  and  performing  angular  scans  across  the  ion  cloud  to  provide 
a  two  dimensional  map  of  ion  column  density. 

2.  Instrumental  Parameters 

The  lar.^e  aperture  (100  mm)  FP1^  achieved  an  overall 
finesse  of  20  and  was  coupled  to  a  servo  controlled  16"  Cassegrain  tele¬ 
scope.  The  telescope  reduced  the  field  of  view  of  the  FP  to  a  point  where 
small  scale  structure  (~  100  meter)  could  be  resolved,  while  retaining  an 
acceptable  light  gathering  power.  Coupled  to  the  telescope  pointing  axis 
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were  an  image  orthicon  television  camera  for  boresight  and  tracking  records 

i- 

and  a  filter  photometer  (monitor)  for  recording  integrated  line  intensities. 
Fields  of  view  were:  for  the  FP-telescope-0. 76  m  rad,  for  the  monitor— 1.66 
m  rad,  and  for  the  TV-  approximately  6°  by  8°.  The  FP  wan  capable  of 
being  pressure  scanned  in  wavelength,  and  achieved  lock  on  the  hyperfine 
line  by  closing  the  inlet  gas  valve  at  a  selected  chamber  pressure.  The 
dual  detecting  head  permitted  observations  of  two  lines  simultaneously, 
normally  X4934  and  X6l42.  The  monitor  could  be  cycled  through  4  filters 
(X4934,  5535,  6111,  and  6l4.2)  and  one  dark  position  for  dark  counting  rate, 
01  could  be  fixed  on  any  one  of  the  above  positions.  It  was  normally  left 
fixed  on  X4934A. 

Normal  operation  of  the  FP  was  in  the  spatial  scan  mode  either  at 
low  speed  to  resolve  fine  structure,  or  at  a  higher  speed  to  record  an 
entire  cloud  profile.  At  release  time  and  at  intermittent  times  throughout 
the  event,  X4934  wavelength  scans  were  taken  to  calibrate  the  optical  depth 
of  the  hyperfine  line.  An  occasional  X5535  wavelength  scan  was  taken  to 
study  the  neutral  cloud  morphology. 

In  the  dual  detection  mode  the  peak  of  the  X6l42  metastable  ion 
line  occurred  at  the  same  FP  chamber  pressure  as  the  peak  of  the  hf  line 
of  X4934.  Thus  a  spatial  scan  in  6l42  light  was  recorded  at  the  same  time 
as  the  hf  scan.  While  the  X6l42  signal  was  not  strong  enough  to  permit 
detection  of  individual  striae,  the  overall  envelope  of  the  cloud  in  X6l42 
will  provide  information  on  the  partition  of  the  ion  population  between 
metastable  and  ground  states. 

Absolute  intensity  calibration  of  the  FP  and  monitor  was  performed 

1  )• 

following  each  event.  Tiie  source  used  was  a  C4-^  radioactive  phosphor  low 
brightness  source  (LBS)  with  calibration  traceable  to  NBS. 
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Real  time  data  was  presented  on  a  Brush  chart  recorder  and  included 
analog  signals  for  chamber  pressure  (wavelength),  FP  intensity  (two  channels), 
monitor  intensity,  differential  azimuth  and  elevation  readouts,  and  several 
channels  of  coded  digital  information,  including  a  one  second  marker 
referenced  to  WWV,  The  same  information  was  also  recorded  in  a  more  accurate 
form  on  a  l/2"  digital  tape  recorder  for  later  computer  retrieval.  The 
television  display,  with  time,  azimuth  and  elevation  superimposed,  was 
recorded  on  a  video  tape  recorder. 

3*  Preliminary  Results 

The  system  was  operational  on  all  six  events  of  the  II  series.  Ho 
information  was  obtained  on  event  Quince  because  of  the  high  ambient  light 
level.  At  this  time  only  the  chart  records  have  been  examined  in  detail; 
however,  it  has  been  determined  that  intelligible  data  do  exist  on  the 
digital  tapes.  Analysis  has  been  concentrated  on  event  Spruce  because  the 
caliber  of  data  is  highest  there  and  because  it  was  generally  the  most 
interesting  event  from  the  standpoint  of  this  system. 

a.  Spruce 

The  most  pertinent  question  that  has  arisen  in  the  II  series 
observations  concerns  the  extent  of  the  optical  thinness  of  the  hf  line  of 
\493^«  On  the  Alaskan  event  Gum  the  lowest  ratio  of  main  to  hyperfine 
components  (main/hf  ratio)  recorded  was  2.1  at  R  +  5  min.  A  first  order 
analysis  (uniform  source  excitation)  related  this  value  to  an  optical  depth 
near  or  below  unity  for  the  hf  line  (see  Section  II).  All  the  II-series 
events,  however,  gave  main/hf  ratios  as  low  as  0.3  and  the  ratio  generally 
remained  at  or  below  1.0  during  most  of  the  visible  cloud  lifetime.  Table  II 
shows  the  values  recorded  during  event  Spruce  along  with  the  spatial 
locations  at  which  the  line  profiles  were  taken.  The  regime  of  ratios  with 


Table  11 


Time 

(after  release) 


hf  half -width  aain/hf 

(®X)  ratio 


location 


1  min. 

100 

0.2 

lb  rain. 

71 

0.32 

19  ain. 

68 

0.33 

23  min. 

30  sec. 

76 

0.22 

23  min. 

40  sec. 

71 

0.55 

24  min. 

20  sec. 

71 

1.1 

25  min. 

70 

0.9 

29  rain. 

20  sec. 

76 

3.4 

30  min. 

76 

1.1 

cloud  center 
solidifying  striation 
aein  body 

striation  on  eastern  edge 
void  on  eastern  edge 
western  edge 
northern  nebulosity 
off  cloud 


on  cloud 
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values  below  unity  can  not  be  treated  by  a  first  order  analysis.  One  of 
the  fundamental  tasks  these  data  call  for  is  a  more  thorough  examination 
of  the  radiation  transport  problem  to  discover  the  significance  of  these 
highly  absorbed  profiles.  One  clue  may  be  that  in  Alaska  the  geometry  was 
such  that  the  X^93^  line  was  observed  to  seme  extent  in  reflection,  while 
in  Florida  the  line  was  seen  much  more  in  transmission. 

Another  measure  of  optical  depth  effects  is  the  amount  of  broadening 
of  a  spectral  feature.  It  is  instructive  to  inspect  the  profiles  from 

this  viewpoint.  Figure  12  shows  two  profiles  taken  at  early  and  middle 
times  in  the  Spruce  cloud’s  visible  lifetime.  Also  shown  as  solid  lines 
are  optically  thin  profiles  including  instrumental  broadening  which  have 
been  adjusted  to  fit  the  hf  line  in  half -width  and  height.  The  temperature 
necessary  to  provide  the  appropriate  half -width  for  the  early  profile  is 
4000 °K.  This  temperature  is  extremely  large,  even  considering  the  evidence 
suggesting  that  the  cloud  remains  at  an  elevated  temperature  (2000°K  for 
Gum)  at  times  as  late  as  one  minute.  In  addition  the  profile  itself 
suggests  absorption  effects,  but  not  by  as  much  as  an  order  of  magnitude. 

If  indeed  there  is  absorption  the  error  induced  is  probably  not  more  than 
a  factor  of  two  even  at  such  early  times.  The  1500°K  optically  thin  profile 
provides  a  surprisingly  good  fit  to  the  lower  experimental  profile.  1500°K 
is  high  but  not  completely  out  of  the  question  for  these  altitudes.  A 
Jacchia-type  estimate  of  exospheric  temperature  at  that  time  should  be  made 
to  determine  just  how  unreasonable  1500 °K  is.  Table  II  shows  that  the  hf 
half -width  does  not  change  appreciably  over  a  wide  range  of  times  and 
main/hf  ratios  in  Spruce,  so  perhaps  what  we  are  seeing  in  the  lower  profile 
is  a  high-temperature,  optically  thin  hf  line. 
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Pending  further  study  we  must  leave  unanswered  the  question  of  the 
optical  depth  of  the  hf  line.  The  remainder  of  this  section  deals  with 
the  spatial  characteristics  of  the  Spruce  cloud  as  viewed  from  site  3a. 

Figure  13  shows  four  successive  spatial  scans  at  medium  speed 
across  a  portion  of  the  Spruce  cloud  covering  times  from  00:ll:58z  to 
00:12:55Z.  The  element  of  emphasis  _n  this  figure  is  the  striation  indi¬ 
cated  by  the  dashed  vertical  line.  It  has  a  half  width  (depending  upon 
choice  of  base  line)  of  less  than  the  resolution  of  the  instrument  or  less 
than  about  150  meters  assuming  l8o  km  slant  range.  For  a  strong,  narrow 
striation  such  as  this,  which  showed  quite  distinctly  on  the  TV  presentation, 
the  modulation  depth. 


%  = 


I  -  I 
peak  valley 

1valley 


(2) 


is  k0%  in  trace  1,  and  65$  after  the  striation  has  expanded  in  trace  4. 
This  series  is  a  particularly  good  example  of  a  well  defined  fine  scale 
striation.  In  most  cases  the  modulation  depth  for  fine  scale  structure 
is  closer  to  10$.  What  appear  visually  as  distinct  striae  are  virtually 
alt  rays  superimposed  on  a  diffuse  background  itfien  they  occur  in  what  was 
once  the  main  body  of  the  ion  cloud. 

Figure  14  is  an  example  of  a  different  class  of  effect  which 
occurred  quite  often  although  not  often  so  dramatically.  The  solid  line 
is  the  hf  peak  intensity  on  a  scan  through  the  late-time  Spruce  cloud, 
and  the  dashed  line  is  taken  from  the  monitor  at  the  same  point.  The 
curves  have  been  normalized  at  the  local  minimum  near  the  center,  and 
they  track  reasonably  well  going  west  from  that  point.  But  the  sharp 


rise  in  the  hf  intensity  going  east  is  not  matched  at  all  by  the  monitor. 

Thus  an  instrument  such  as  a  camera  or  the  monitor  which  records  integrated 
line  intensity  would  give  an  erroneous  value  for  the  ion  column  density. 

The  same  effect  quite  often  shows  up  when  scanning  across  a  localized 
enhancement  such  as  a  striation.  The  modulation  depth  is  greater  in  the 
hf  line  than  in  the  total  line  as  recorded  by  the  monitor. 

One  of  the  stated  objectives  of  this  system  was  to  track  the  flashing 
light  on  the  rocket  payload  vehicles  (probe  and  beacon)  on  the  TV  to  obtain 
a  point  by  point  comparison  of  RF  electron  densities  with  optical  ion 
densities.  Only  two  vehicles  have  thus  far  been  identified  on  video  tape, 
both  on  event  Spruce.  Table  III  gives  times  during  which  these  vehicles 
were  observed  and  the  operational  mode  of  the  FP  surrounding  those  times. 
Exact  vehicle  positions  will  be  taken  from  the  video  tape. 

During  most  of  the  beacon  occultation  as  viewed  from  site  3  the 
FP  was  performing  X4934  wavelength  scans  on  a  developing  striation  near 
the  trailing  edge  of  the  ion  cloud.  One  of  the  scans  is  shown  in  the 
lower  trace  of  Figure  12.  At  the  time  of  the  last  observed  flash  the 
beacon  was  traversing  that  same  striation  at  a  point  about  2°  south  of 
where  we  had  taken  the  wavelength  scan.  It  was  some  3 A  °  into  the  cloud 
on  a  track  running  at  approximately  65°  to  the  striation  axis. 

Figure  15  shows  the  hf  line  traces  of  the  three  spatial  scans 
surrounding  the  occultation.  The  earliest  of  the  three  (top)  shows  the 
abrupt  change  in  slope  associated  with  the  hard  edge,  and  the  bridge 
trailing  to  the  left  in  the  figure.  The  middle  trace  shows  sane  motion 
of  the  cloud  to  the  left  (eastward),  and  very  little  structure.  Fine 
structure  could  be  present  but  the  chart  records  do  not  have  enough 


Table  III 


Event  Spruce— Times  surrounding  vehicle  observations  on  TV. 


Time  (UT) 


Event 


23:57:51 

58:18 

58:28 

58:45 

59:02 

59:03 

59:23 

59:30 

59:38 

00:00:07 

00:35 

01:02 


Begin  spatial  scan* 
end  spatial  scan 

begin  spatial  scan-through  bridge  and  across  long 
axis  of  ion  cloud 
acquire  probe 
lose  probe 
end  spatial  scan 
reacquire  probe 
last  observed  probe  flash 
begin  spatial  scan 
end  spatial  scan 
begin  spatial  scan 
end  spatial  scan 


00:03:55 
04: 08 
04:19 
04:46 
04:55 
05:39 
05:42 
05:50 
06:08.3 
06:21 
06:22 
06:25 
06:46 
07:07 
07:46 
08:00 
08:23 


begin  spatial  scan 
end  spatial  scan 
begin  spatial  scan 
end  spatial  scan 
begin  spatial  scan 
end  spatial  scan 
acquire  beacon 

begin  X493^  wavelength  scan  on  solidifying  striation 

beacon  passes  visible  hard  edge 

end  X4934  wavelength  scan 

last  observed  beacon  flash 

begin  X5535  wavelength  scan  on  neutral  cloud 

end  X5535  wavelength  scan 

begin  spatial  scan 

end  spatial  scan 

begin  spatial  scan 

end  spatial  scan 


*unless  otherwise  noted,  the  spatial  scans  are  approximately  parallel 
to  the  long  axis  of  the  ion  cloud  through  its  center  of  gravity. 
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definition  to  be  positive.  The  vertical  bar  shows  the  farthest  excursion 
from  the  smoothed  curve.  The  bottom  trace  (post-occultation)  shows 
structure  which  is  definite  but  not  extremely  pronounced. 

At  the  present  time  absolute  ion  densities  have  not  been  obtained. 

An  absolute  calibration  using  the  phosphor  LBS  was  taken  following  the  event 
but  needs  to  be  retrieved  from  the  digital  tape.  Also  an  estimate  of  optical 
depth  effects  should  be  made  before  the  densities  can  be  relied  upon. 

b.  Other  Events 

Data  reduction  ha*  been  concentrated  on  event  Spruce,  but 
some  observations  have  been  made  about  the  general  character  of  other 
releases. 

Small  scale  (less  than  the  FP  resolution  element)  structure  was 
observed  on  all  releases  except  Event  Olive.  Striation  was  seen  on  Olive 
but  appeared  to  be  the  preliminary  large  scale  (~  1  km)  modulation  that 
characterized  the  other  clouds  before  fine  scale  structure  appeared.  The 
early  \h93^  profiles  on  Olive  showed  line  reversal  at  the  peak  of  the  hf 
line  and  should  provide  a  good  test  of  the  transport  theory. 

On  events  Nutmeg;  Plum  and  Spruce  the  large  scale  striation 
structure  eventually  divided  the  smooth  envelope  of  the  ion  cloud  into 
discrete  groups  of  striations,  giving  a  very  irregular  profile,  as  shown  in 
Figuies  13  and  14.  This  happened  before  the  small  scale  structure  became 
apparent.  On  Redwood  such  was  not  the  case;  the  small  scale  structure 
evolved  as  modulation  on  the  smooth  ion  cloud  envelope,  whose  general 
character  remained  unchanged  throughout  its  visible  history.  Modulation 
depths  on  the  main  body  were  typically  10^-20$.  Event  Olive  gave  indica¬ 
tions  that, its  evolution  was  following  a  Spruce-like  course,  but  on  a 
much  longer  time  scale. 
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Release  times  of  the  main  cloud  and  two  of  the  three  puffs  have 
been  taken  from  video  tape  on  Olive.  The  main  cloud  appeared  at  23:53:57»3Z 
and  was  followed  by  puffs  at  23:54:06. 2Z  and  23 .*54: 15. 9Z,  which  both  merged 
into  the  main  cloud. 

The  third  puff  which  evolved  into  a  separate  cloud  showed  distinct 
small  scale  structure  at  9  minutes  after  the  main  cloud  release.  This 
small  scale  structure  was  gone  at  R  +  18  and  only  the  large  scale  modulation 
remained. 

The  few  X5535  wavelength  scans  taken  during  events  Olive  and 
Spruce  showed  strong  reversal  immediately  after  release,  but  relaxed  to  a 
more  thermal  line  shape  at  middle  times  in  the  cloud  history.  The  early 
time  line  shapes  will  provide  an  excellent  test  for  the  hypothesis  developed 
from  Alaskan  event  Gum  that  the  cloud  remains  at  an  elevated  temperature  for 
some  time  after  release. 

4.  Data  Reduction 

Currently,  efforts  have  been  underway  to  retrieve  the  data  from  the 
digital  tapes.  The  form  the  data  bank  will  take  is  planned  to  be  a  listing 
for  each  spatial  scan  giving  azimuth,  elevation,  differential  angle,  FP 
intensity,  and  monitor  intensity  versus  time.  Requests  for  specific  time 
scans  can  then  be  satisfied  with  a  copy  of  the  listing,  or  a  computer 
drawn  plot  if  possible.  The  wavelength  scans  will  be  presented  as  FP 
intensity  versus  pressure  (wavelength). 

It  is  planned  to  correlate  certain  spatial  scans,  such  as  in 
Figure  l3swith  photographs  taken  up  the  field  lines  to  arrive  at  a  three 
dimensional  ion  density  model. 


APPENDIX  A 


RESONANCE  RADIATION  TRANSPORT 

For  interpretation  of  the  bariun  profiles  recorded  in  Alaska, 
it  has  been  necessary  to  develop  a  computer  code  to  generate  line  profiles 
from  a  model  of  the  bariun  cloud.  For  the  theory  of  resonance  radiation 

transport  in  a  spherically  symmetric  geometry  we  nave  drawn  extensively 

5 

on  the  paper  by  Thom«6.  This  Appendix  details  the  application  of 
Thomas'  theory  to  our  problem. 

As  mentioned  in  the  text  our  approach  to  the  bariun  neutral  cloud 
is  to  treat  X5535  transport  as  a  two  level  process,  with  the  upper  level 
of  the  X5535  line  excited  only  by  capture  of  a  X5535  photon,  and  decaying 
only  through  the  same  channel.  Thus  the  neutral  metastable  states  are 
ignored,  as  "ell  as  cascade  to  and  excitation  fran  the  upper  level.  This 
approximation  is  probably  quite  good  in  late  times  when  the  greatest 
optical  thickness  of  the  cloud  is  about  3>  "but  becomes  more  questionable 
in  early  times  when  the  optical  thickness  is  much  higher.  Further  work 
is  required  to  determine  to  what  extent  competing  processes  contaminate 
the  two  level  system. 

An  approximation  inherent  in  Thomas'  treatment  is  that  of  complete 
frequency  redistribution  (CFR),  which  is  that  the  emission  coefficient 
and  absorption  coefficient  are  equal.  A  further  approximation  of  Thomas 
is  that  both  emission  and  absorption  coefficients  are  Gaussian.  Both  of 
these  approximations  are  tied  together  since  CFR  is  commonly  assuned  to 
be  valid  (see  Thomas)  in  the  dopplercore  of  the  line  (out  to  3  doppler 
widths  from  the  line  center).  It  is  near  this  same  point  that  the 
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magnitude  of  the  natural  line  profile  of  the  X5535  lint  (natural  half¬ 
width  0.3  mK,  doppler  half-width  33  mK)  hecomeB  comparable  to  the 
magnitude  of  the  doppler  profile.  So  both  approximations  may  be 
accepted  after  one  minute  into  the  cloud  history  when  the  line  width 
is  2 4  doppler  widths  or  less.  Earlier  than  that,  both  approximations 
are  suspect,  but  extremely  laborious  to  investigate.  In  particular  it 
would  be  desirable  to  see  if  a  Voigt  (doppler  plus  natural)  emission 
profile  (used  in  place  of  a(x)  in  Eq.  (A'.-lO))  would  supply  enough  signal 
in  the  wings  of  the  line  to  affect  the  conclusions  derived  from  Figs.-  3  &  4 
i.e,  that  a  high  temperature  is  required  to  reproduce  the  observed 
profile  at  R  +  45  sec. 

Thomas  has  expressed  the  equation  of  radiative  transfer  in  the 
form  of  an  integral  equation  for  the  excited  state  density  (Eq.  (35)  in 
Thomas’  paper): 

N(r)  »  N(l)(?)  +  k0p(r)  f  ^  ds  G(r,?>(r#)  ,  (A.l) 

where  N^(r)  is  the  density  of  atoms  excited  by  the  primary  solar 
radiation: 

(l)  k 

N  (r)  =  (J0/hv0Y)k0p(r)TT  dvD  T(r,«,JTg)  ,  (A.  2) 

with  JQ  the  solar  flux  per  unit  frequency  (assuned  constant  over  the  line)j 
h  -is  Planck’s  constant,  vQ  is  the  center  frequency  of  the  line,  y  is  the 
damping  width,  kQ  is  the  peak  absorption  cross  section  p(r)  is  the 
spherically  symmetric  density  distribution  of  ground  state  atoms,  dv^  is 
the  thermal  doppler  width  and  T  accounts  for  the  attenuation  of  the 
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primary  solar  flux  before  reaching  F>  T(r,r/)  is  the  probability  that 
a  photon,  emitted  at  the  point  F,  will  travel  to  the  point  r*  without 
being  absorbed. 

T(r,r')  =  rr"?  J  a(x)exp(-a(x)T(r,r')dx  (A. 3) 

«CP 

where  x  is  a  dimensionless  frequency  variable 

X  =  (v-vq)/Avd  (i.4) 

and 

or(x')  =  exp(-x2)  (A.  5) 

is  the  (normalized)  absorption  and  emission  coefficient  (equal  because 
of  CFR),  and 

T(F,r')  =  kQ  JL,  p(r")ds"  (A.6) 


is  the  optical  depth  between  r  and  r',  with  the  integration  being 
carried  out  in  a  straight  line.  T(r,e°>n_)  then  is  the  probability  that 
a  primary  solar  photon  incident  upon  the  mediua  from  a  direction  will 
reach  F  without  being  absorbed. 

The  second  tern  in  Eq.  (A.l)  represents  radiation  emitted  at  F# 
and  absorbed  at  r.  The  integration  ds~  is  carried  along  a  line  in 
direction  R,  with  R  integrated  over  all  possible  directions.  Gfor7)  is 
the  (unnormalized)  probability  that  a  photon  emitted  at  F/  will  be 


absorbed  at  r. 


k2 


G(r,r')  =  tr**f  [«(x)]2  exp(-ar(x)T(r,r/)]<3x  .  (A-7) 

The  expressions  used  for  calculating  T  and  G  are 

T(t)  *  S*(-T)n/n.,(n+l)^  0  St<« 

°*  T-1(n  lu  t)~*  3.0  <T<«  (A.8) 

and 

G(T)  «  E"(-T)n/n.f(n+2)^  0  <f  <» 

•*  (1  +  \tn  t)t“2  (tt  /n  t)"1  8.0  s  t  <  •  ,  (A.9) 

where  the  approximations  were  used  when  the  series  became  unstable  due  to 
computer  round-off  error. 

Once  N(r)  has  been  determined,  a  line  profile  may  be  obtained  by 
integrating  along  the  line  of  observation 

os 

.  l(x,ZT  )  =  constant  x  f  cr(x)  exp[-»(x)T(r/,®,iT  )N(r/)ds  (AvlO) 

where  the  magnitude  of  the  constant  is  not  important  since  we  are 
concerned  only  with  shape. 

Iterative  techniques  are  used  to  solve  (A*l)  for  N(r).  Two 
problems  arise  in  this  connection.  The  first  is  that  convergence  using 
N^(r)  as  a  trial  solution  is  a  prohibitively  long  process,  since 
N^(r)  is  generally  very  far  from  the  solution  N(r) .  This  problem  is 
dealt  with  by  using  an.  approximation  .(c.f.  Thomas)  to  obtain  a  more 


*8 

accurate  trial  solution.  If  N(r/)  in  Eq.  (A.l)  is  assuned  constant  over 
the  region  where  G(r,r')  is  appreciable  and  taken  out  of  the  integration 
sign,  and  the  equation  simplified  by  means  of  an  analytic  relation  between 

T  and  G,  it  .becomes 


M(r)  =  M^(r)/B(r) 


(A.  11) 


where 


E(r)  =  J'(dfi/4rr)T(r,»)  =  1  -  kQ  J  dfi/hn  J  d?  p(r')G(r,7') 

(A.  12) 

is  simply  the  probability  that  a  photon  emitted  at  point  r  escapes  the 

mediun.  completely  without  being  reabsorbed*  N^(r)/E(r)  is  then  used 

as  a  first  approximation  for  the  iterative  solution  of  (B.l).  This 

approximation  along  with  the  preceding  transport  theory  is  discussed  by 

5 

Thomas . 

Our  final  solution  does  not  use  the  iterative  technique  inherent 
in  (A.l).  Instead  we  have  used  Eq.  (A. 12)  and  the  Born-Oppenheimer 
principle  of  expansion  of  an  integral  equation  to  rewrite  (A.l)  in  the 
series  form: 

H(r)  =  M^(r)/E(r)  +  kQ  J'  (dfi/Jm)  J'  ds[p(r)H^(r')/E(r#) 

-  p(r')N(l)(F)/E(r)]  x  G(r,r')  +  1^(7)  +  N4(7)  + .  (A.13) 

where  N^(r)  and  successive  terms  are  related  to  the  immediately  preceding 
term  by: 


- 


t  a 


!  'S 


1 

1 


is 
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Nk(7)  -  k0  p(r)  j*  dfi/4tT  J  ds  N(k  l)(?/)G(F,?')  .  (A.  14) 

The  main  attribute  of  Eq.  (A. 13)  is  that  its  leading  term  is  close  to  the 
final  solution  and  it  was  felt  that  nuuerical  errors  in  the  calculation 
of  G  would  not  propagate  to  the  same  extent  as  in  Eq.  (A.l). 

Convergence  requires  roughly  .  .half .  as  many  terms  as  the 
mediun  is  optically  thick,  if  the  series  is  terminated  when  the  last  term 
adds  less  than  0.1$  of  the  sua  of  the  preceding  terms  at  all  mesh  points. 

A  spherical  coordinate  system  was  chosen  with  the  origin  at  the 
cloud  center  and  the  polar  (0=0)  axis  along  the  sun-cloud  line  ( UB ),  as 
shown  in  Figure  9.  Since  the  solar  flux  is  symmet. ' 2  about  the  polar 
axis,  the  excited  state  distribution  will  also  have  cylindrical  symmetry 
about  that  axis,  N(r)  =  K(r, 0).  Therefore  we  need  only  find  N(r)  for 
0Sr<“,  0  s  0  S  tr.  The  radial  mesh  interval  chosen  was  Ar  =  g/4  where 
g  is  the  folding  distance  of  the  spherically  symmetric  Gaussian  ground 
state  distribution 

p(r)  =  PQ  exp(-r2/g2)  ,  (A-15) 

and  the  cloud  was  terminated  at  n  =  3g«  The  8-interval  was  chosen  to  be 
tt/8,  as  was  the  interval.  For  mmerical  solution  the  value  of 
N(r,  0)  was  determined  for  each  r  and  0  mesh  point.  Equation  (A.l4) 
involves  a  three  dimensional  integration,  but  ^ (r/)  is  independent 
of  cp#,  while  G(r,r/)  still  is  tp/  dependent.  So  the  cp/  integration 
(carried  around  the  torus  indicated  in  Figure  9)  may  be  incorporated 
into  G  and  (A.l4)  rewritten 


_  IJ 
.2  .  G 

I  ,J  T'J' 


iV 

w 

(k-i) 


(A.  16) 


with  I  and  J  representing  the  radial  and  angular  variables.  To  obtain  G 
each  volume  element  was  further  divided  into  8  subelements.  But  it 

IJ 

was  still  necessary  to  find  the  self-contribution  of  ^ach  elooieixi  G^j, 
not  frcm  direct  integration,  but  by  using  Ea.  (A-12)  to  write 


2 


(A*  IT) 


and,  in  a  sense,  normalizing  the  self-contribution  to  account  for 

numerical  mistakes  in  the  other  elements  and  to  conserve  photons. 

An  example  of  the  excited  atom,  distribution  obtained  is  shown 

in  Figure  10.  Even  though  the  full  series  solution  N(r)  differs  frcm 

'1)  . 

the  first  order  approximation  N'  /E  along  the  sun-cloud  line,  they 
track  almost  perfectly  through  the  center  at  right  angles  to  that  line. 
Since  our  observations  were  made  nearly  at  right  angles  to  the  sun-cloud 
line  we  were  able  to  save  considerable  computing  time  by  using  the 
IT^/E  approximation  to  generate  Figure  3*  It  was  necessary  to  find 
the  final  solution  to  be  sure  trust  the  approximation  was  valid. 

For  handling  optical  depth  effects  on  she  X^93^  ion  line  we 
desired  only  to  find  a  first  approximation,  since  the  cross  excitation 
between  ground  and  metastable  levels  is  much  more  important  than  in  the 
neutr\l  case,  ana  since  the  rear  geometry  is  completely  intractable. 

The  approximation  we  have  chosen  is  to  assume  that  radiation  transport 
will  supply  an  excited  ion  distribution  proportional  to  the  ground  state 


46 


/ 


distribution,  and  the  only  optical  depth  effects  ccme  from  reabsorption 
of  emitted  radiation  before  it  reaches  the  detector.  One  justification 
for  this  approximation  is  that  in  applying  this  crude  model  to  the 
neutral  cloud  it  supplied  a  half-width  that  agreed  with  the  full  transport 
solution  to  within  15 over  a  range  of  optical  depths  from  0.01  to  1000. 

} 

An  increase  in  the  optical  depth  of  the  median,  is  manifested  by  a  ! 

decrease  in  the  ratio  of  the  main  to  hyperfine  component  from  its 

j 

optically  thin  value  of  lk.7  (see  Figure  7).  We  have  attempted  to 
relate  this  ratio  to  the  optical  depth  (ion  .line  of  sight  density). 

To  implement  this  approximation  Eq.  (A.IO)  is  rewritten  in 
terms  of  optical  depth  instead  of  distance,  dr  =  k0p(r)ds; 

x 

)  -  const  .  j*  a(x)  exp[-ff(x)T,]N(T#)dT/  (A. id); 

0  b 

where  N(t/)  =  N(F#)/(kop(r#)).  At  line  center  xeO  and  using 
N(t#)  =  constant  (N(t)  «  p(t)) 

I  =  const  x  N(t)[1  -  exp(-t)]  .  (X.19)  | 

! 

5 

1 

| 

Now,  since  both  main  (m)  and  hyperfine  (hf)  level  populations  are  assuned  j 

1 

proportional  to  the  ground  state  population,  Nm(T)  *=  14.  7  Nhf(r)  and  j 

Tffi  *=  14,7  we  may  take  the  ratio  of  main  to  hyperfine  peaks;  j 

i 

i 

i 

Iom/Iohf  c  Cl  -  exp(-T))/[l4.7(l  -  exp(-r/l4.7)) j  (A. 20)  j 

t 

where  now  T  refers  to  the  optical  depth  of  ions  responsible  for  absorbing 


the  main  component,  i.e.  82$  of  the  total  in  the  ground  state.  The 
predicted  main  to  hyperfine  ratio  (Eq..  (A.  20))  has  been  plotted  in 
Figure  11  against  T  which  is  defined  as  the  optical  depth  of  al3  ions 
in  the  ground  state,  i.e.  T  «  t/,82.  Figure'  U  has  then  been  used  for 
our  reduction  of  the  ion  line  data. 


In  the  late  stages  of  this  research  an  investigation  was  made 
into  using  emission  and  absorption  profiles  which  are  more  general  than 
the  Gaussian  profiles  initially  used  in  the  neutral  barium  line  shape 
calculations.  It  was  not  feasible  to  include  natural  line  broadening  in 
the  transport  section  of  the  calculations  (determining  the  spatial 
distribution  of  atoms  in  the  upper  level),  although  it  is  planned  to  do 
thi6  in  the  future.  However,  the  appropriate  Voigt  (natural  plus 
doppler  broadening)  profile  has  been  used  for  the  emission  and  absorption 


profiles  which  enter  the  line-of-sight  integration  for  the  recorded 
profile.  That  is,  or(x)  in  Eq.  (A.  10)  has  been  replaced  by  the  Voigt 


and  "a"  is  the  ratio  of  natural  to  doppler  line  widths.  Figs.  3> 

and  5  have  been  produced  using  the  results  of  these  hybrid  calculations. 
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Figure  Captions 

(1)  Sample  of  real  time  chart  records  obtained  during  barium  releases. 
These  scans  show  1.493^  and  15535  lines  from  the  early  moments  of 
event  C-um.  The  information  in  the  indicated  channels  is  as  follows: 
Channel  information 

1,  6  Interference  filter  code  (not  operational  here) 

2  Gamma-scope  memory,  showing  buildup  of  lines 

3,  8  10—  channel  marker 

4  Pressure  transducer  with  filter  change  and  start  signals 
superimposed 

5,  10  Serial  time  code  (GMT,  erratic  operation  here) 

7  Monitor  scalar  (photometer) 

9  FP  scaler  (single  pass  line  profiles) 

(2)  Partial  term  diagrams  for  neutral  and  ionized  barium,  with  the 
transitions  we  have  studied  listed  by  wavelength  in  A.  Also  shown 
are  the  metastable  levels  and  some  levels  which  link  the  meta¬ 
stables  together  and  to  the  ground  level. 

(3)  Observed  15535  contrast  finesse  for  barium  release  GUM,  and 
curves  showing  predicted  variation  of  contrast  finesse  for  a 
cloud  of  constant  inventory  evolving  by  diffusion,  labelled  by 
the  total  Ba  in  the  cloud. 

(4)  Experimental  15535  profile  at  R  +  45  seconds  for  event  GUM,  and 
predicted  profiles  for  different  comoi nations  of  barium  temperature 
and  inventory.  Only  half  of  the  (symmetric)  profiles  are  shown 
and  the  experimental  profile  has  been  folded  about  its  midpoint 

to  enhance  the  comparison. 
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(5)  Time  variation  of  X5535  in  event  GUM.  Vertical  axis  is  relative 
intensity.  Model  curves  have  been  normalized  to  experimental 
curves  at  R  +  2  min. 

(6)  Time  variation  of  X5826  in  event  GUM.  Vertical  axis  is  relative 
intensity.  Smooth  curves  are  visual  fits  to  data  points.  Time 
is  GMT(Z). 

(7)  XU934  profile  observed  at  R  +  26  minutes  in  event  GUM,  and 
calculated  profile  for  no  absorption.  The  arrows  indicate  centers 
of  gravity  of  closely  grouped  lines  and  their  relative  intensities. 
The  peak  at  ~75  on  the  frequency  scale  is  the  "hyperfine  line". 

(8)  Time  variation  of  Xl93**  in  event  GUM.  Vertical  axis  is  relative 
intensity,  time  is  GMT(z).  Subsidiary  component  is  the  "hyperfine" 
line. 

(9)  Resonance  radiation  transport  geometry.  Since  the  exciting  solar 
flux  (incident  along  Qs)  is  symmetric  about  the  z-axis,  as  is  the 
ground  state  Ba  density,  then  the  excited  state  density  will  be 
also.  Therefore,  a  solution  need  only  be  found  in  a  semicircular 
plane,  shown  hatched  in  the  plane  of  the  figure.  The  three- 
dimensional  path  which  must  be  employed  to  find  the  transfer 
matrix  G(r,r/)  is  indicated  by  an  ellipse.  The  line  of  sight 
(along  Qc)  integration  for  the  line  profile  is  carried  from  one 
edge  of  the  cloud  to  the  other,  through  the  origin  and  perpendicu¬ 
lar  to  the  z-axis  (fl0  lns).' 

(10)  Excited  state  densities  along  ns(||).  and  along  no(]J  (see 

Figure  9)*  is  the  density  excited  by  primary  solar  flux 

(absorbed  but  not  scattered),  is  'the  Sobolev  (see  text) 

approximation  to  the  final  solution,  and  N  is  the  final  solution. 
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(11)  First  approximation  to  optical  depth  effects  on  X4934.  The 
curve  gives  the  value  of  line  of  sight  ion  ground  state  density 
which  would  produce  a  given  main/hyperfine  ratio  in  X4934,  under 
the  uniform  excitation  approximation  (N(r)  «  p(r)). 

O 

(12)  Solid  points:  X4934  A  line  profiles  obtained  during  event  §pruce 
(R  =  23:52 :04z).  Solid  curve:  computer-generated  optically  thin 
profiles  adjusted  to  fit  the  X4934  hf  line  in  height  and  half¬ 
width.  The  unabsorbed  peak  intensity  is  100  units. 

(13)  Successive  spatial  scans  in  hf  light  over  a  portion  of  the  Spruce 
ion  cloud  showing  a  fine- scale  striation.  Scans  from  left  to  right. 
(1)  00:11:59  to  12:65Z,  (2)  12:13  to  12:19,  (3)  12:32  to  12:38, 

(4)  12:49  to  12:55. 

(14)  Spatial  scan  in  hf  light  and  in  total  X4934  light  of  Spruce 
leading  edge  showing  shape  differences  due  to  optical  depth  effects. 
Scans  from  left  to  right  from  00:19:37  to  19:50Z. 

(15^  Three  spatial  scans  nearest  in  time  to  beacon  occultation  on  event 
Spruce.  Scan  was  performed  from  left  to  right. 
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